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Nanodots with magnetic vortices have many potential apphcations, such as magnetic memories 
(VRAMs) and spin transfer nano-oscillators (STNOs). Adding a perpendicular anisotropy term to 
the magnetic energy of the nanodot it becomes possible to tune the vortex core properties. This 
can be obtained, e.g., in Co nanodots by varying the thickness of the Co layer in a Co/Ft stack. 
Here we discuss the spin configuration of circular and elliptical nanodots for different perpendicular 
anisotropics; we show for nanodisks that micromagnetic simulations and analytical results agree. 
Increasing the perpendicular anisotropy, the vortex core radii increase, the phase diagrams are 
modified and new configurations appear; the knowledge of these phase diagrams is relevant for the 
choice of optimum nanodot dimensions for applications. MFM measurements on Co/Pt multilayers 
confirm the trend of the vortex core diameters with varying Co layer thicknesses. 
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I. INTRODUCTION 

Nanoscopic and mesoscopic magnetic structures have 
attracted the interest of many workers in recent years in 
view of their very interesting physical properties and for 
their potential applications. Quasi- twodimensional mag- 
netic nanodots made of soft magnetic materials, such as 
permalloy, may present for their lowest energy state, sev- 
eral magnetic configurations: (i) quasi-uniform in-plane 
state, (ii) quasi-uniform out-of-plane state, (iii) magnetic 
vortices or swirls^. 

Magnetic vortices are structures where the magnetic 
moments are tangential to concentric circles. The center 
of the vortex has a singularity (the vortex core) where 
the magnetization points out of the plane, with a ra- 
dius, in the thin dot limit^, of th e order of the exchange 
length of the material l^x = \/274//io^s , where A is 
the exchange stiffness constant and Mg is the satura- 
tion magnetization^^^. With the parameters used in the 
present work, for permalloy /ex = 5.3 nm and for cobalt, 
Iqx = 4.93 nm. 

Magnetic vortices have been observed by many exper- 
imental techniques, such as magnetic force microscopy^, 
X-ray microscopy^, or inferred from hysteresis curves^; 
they also result from theoretical modeling^^^. 

The proposed applications of magnetic nanodots 
include their use in patterned magnetic recording 
media^^, as elements in magneto-resistance random 
access memories (MRAM's}^^!^, spin transfer nano- 
oscillators (STN0's)i2i"— and nanoscopic agents for can- 
cer treatment^. 

The different magnetic configurations, and conse- 
quently the large variation in the magnetic properties ob- 
served in nanodots as a function of dimensions, underline 
the interest in the study of diagrams (phase diagrams) 



mapping the parameter space where a given magnetic 
behavior is to be expected. 

Experimental studies have been used to obtain the 
phase diagram for magnetic disks: Ross et ali^ de- 
rived the phases from hysteresis curves, Chung et al^ 
from SEMPA measurements. Metlov and Guslienko2^ 
obtained a phase diagram with regions of in-plane magne- 
tization, perpendicular magnetization and vortex struc- 
ture; the equilibrium magnetic configuration obtained 
by micromagnetic simulation showed general agreement 
with this diagran>^i2^. Another simulation study, this 
time using a scaling approach, was made for circular and 
elliptical nanodots^^. They have found in the phase di- 
agram for ellipses a double vortex arrangement for dots 
with semi- axis a larger than 150 nm. Their simulations 
were made for core- free ellipses, a choice that might dis- 
place the phase boundaries by a significant 35%. 

We have recently shown theoretically and 
experimentally^^, that using Co/Pt multilayers it 
is possible to tailor the vortex core diameter by playing 
with the perpendicular anisotropy originated at the 
Co-Pt interface. When one increases the perpendicular 
anisotropy acting on a magnetic nanodot, e.g., reduc- 
ing the Co layer thickness, the vortex core diameter 
increases, and eventually another vortex state appears, 
which is characterized by an out-of-plane magnetization 
component at the dot rim. The increase in perpendicular 
anisotropy has an effect that is somehow equivalent to 
a reduction in the relative importance of the in-plane 
shape anisotropy. 

The main goal of the present work is to study how the 
phase diagram of magnetic nanodots is modified by the 
presence of perpendicular anisotropy (Kz). However, we 
have first obtained the phase diagram for nanodots with 
Kz = 0, with circular and elliptical shapes. This has 
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been done for two main reasons, first to illustrate and 
validate our methodology, which will be used in sequence 
in this paper, and to verify the effect of the magnetostatic 
energy (responsible for shape anisotropy) on the ground 
state magnetic arrangement. 

We have obtained phase diagrams by micromagnetic 
simulation and analytically, and they are in agreement. 
We have also indicated experimentally, using Magnetic 
Force Microscopy (MFM), that the vortex core diameter 
shows the same trend as predicted by the theory. The 
paper is organized as follows: in section II we discuss 
magnetic configurations of the disks, which include IIA) 
micromagnetic simulations, IIB) analytical description, 
IIC) experimental results. In Section III we present the 
results for the ellipses: IIIA) micromagnetic simulations. 
In Section IV we present a brief discussion, a summary 
of the main results with the conclusions. 



II. RESULTS FOR DISKS 
A. Disks: Micromagnetic Simulations 

For the simulations we used the OOMMF package (a 
free software available from NIST^^), using the parame- 
ters for bulk permalloy, to allow a comparison with the 
literature (exchange stiffness constant A = 1.3 x 10~^^ 
J/m, saturation magnetization Mg = 860 x 10^ A/m). 
We have neglected the in-plane anisotropy; however, we 
have also simulated magnetic systems exhibiting a per- 
pendicular anisotropy. An application of such simula- 
tions is the description of the behavior of the Co/Pt mul- 
tilayers. To make the present results of more general use, 
they have been given in terms of normalized parameters, 
using the exchange length Iq^- 

For some dimensions of the nanodots, the simulation 
may converge to a configuration that does not correspond 
to the absolute energy minimum. Therefore, the simula- 
tions made with parameters near the boundary regions 
between different configurations had to be made by ini- 
tially imposing different magnetic configurations, and af- 
ter the convergence, comparing the resulting energies to 
determine the spin arrangement corresponding to the ab- 
solute energy minimum. The boundary lines between the 
different phases were obtained from the intersection of 
the curves of total energy for the different states^^. 

Effects of discretization are inherent in the methodol- 
ogy used here (see^^). For this reason we have studied 
the effect of the cell size (for most of our simulations 
5x5x5 nm^). We have found that the position of the 
boundaries of our phase diagrams change very little for 
different cell sizes; these effects are even less important 
for the configuration with perpendicular magnetization. 

The phase diagram for magnetic nanodisks obtained 
from the computed energies for the different magne- 
tization configurations and for different perpendicular 
anisotropy values (Kz) is shown in Fig. [H For Kz = the 
phase diagram agrees with those of references?- J^i^5^i24. 




FIG. 1: (Color online) Phase diagrams for circular nanodots 
as a function of reduced height h/lex and reduced radius 
r/lex, drawn from the minimum energy computed with mi- 
cromagnetic simulation for different values of the perpendic- 
ular anisotropy: = 0, 100 x 10^ J/m^ 200 x 10^ J/m^ 
360 X 10^ J/m^. The different magnetic configurations are 
labeled in the inset showing the Kz = phase diagram: a) 
single domain parallel to the plane, b) single domain perpen- 
dicular to the plane, c) vortex, d) configuration given by Fig. 
[2]i, and e) vortex with perpendicular component (Fig. [2^). 



It shows three regions, depending on the aspect ratio of 
the disks. The corresponding magnetic configurations are 
shown in Fig. [21 For very thin disks, for a wide range 
of disk radii, the shape anisotropy favors a quasi-uniform 
in-plane state; on the other hand, for thicker disks and 
approximately r < 4 /gx, a quasi-uniform out-of-plane 
state is observed, that is easy to understand, since in 
this region nanodots cannot be taken as approximately 
twodimensional disks. Finally, for r > 4 l^^ and h > 41^^ 
we observe different magnetic vortex states as the ground 
states. One should note that, in the vortex state region in 
the graph, the magnetization shows an increasing out-of- 
plane component as the disk thickness increases; this can 
be seen in Fig. [2^. Note that the vortex core diameter 
varies along the length of the cylinder, reaching a maxi- 
mum at half the height. Also, we observe the formation 
of a "mixed" state (e in Fig. [1]), where the magnetiza- 
tion shows vortex domains at the dot ends, and out of 
plane magnetization near half height. An axial section of 
a nanodisk with larger h (Fig. [2]i) shows that the vortex 
acquires a perpendicular magnetization component. 

When we include a perpendicular anisotropy term, 
the phase diagram of the disks is modified, as can be 
seen in Fig. [TJ As expected, the region correspond- 
ing to magnetization perpendicular to the plane (region 
b, shown in the inset) is increased as is increased, 
displacing to larger radii the boundary line between the 
quasi-uniform perpendicular moment state and the vor- 
tex state. Furthermore, the region for in-plane magneti- 
zation (region a) is reduced. In the case of the simulation 
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FIG. 2: (Color online) Magnetization configurations for nan- 
odisks: a) quasi-uniform in-plane magnetization; b) quasi- 
uniform perpendicular magnetization; c) magnetic vortex; d) 
disk with parameters in the region above the red continu- 
ous line (or the blue dotted line) in Fig. [3] (d=400 nm and 
anisotropy Kz = 375 x 10^ J/m^); e) lateral view showing 
the longitudinal section of an elongated nanodisk as found in 
the region of the phase diagram where the vortex acquires a 
perpendicular magnetization component (see region e in Fig. 
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FIG. 3: (Color online) Phase diagram for a 10 nm-thick 
magnetic disk, as a function of the perpendicular anisotropy 
Kz and reduced radius r/lex for micromagnetic simulation 
(red continuous line) and analytical computation (blue dotted 
line). The spin configurations are: a) single domain parallel to 
the plane, b) single domain perpendicular to the plane, c) vor- 
tex, d) configuration given by Fig. [2]i. The lines on the right- 
hand side (red continuous=simulation, blue dotted=analytic) 
limit the region above which the spin structure is given by d. 



with the highest perpendicular anisotropy shown in the 
figure {Kz = 360 x 10^ J/m^), the in-plane magnetiza- 
tion region is limited to a narrow range between 7 and 
12 /ecc, for very thin disks. The tendency of an out-of- 
plane magnetization in the vortex region is not observed 
in the cases of nonzero perpendicular anisotropy. For 
higher anisotropics another more complex configuration 
appears with an out-of-plane magnetization at the disk 
rim (region d), as observed experimentally in lef^^. 

Increasing the perpendicular anisotropy increases the 
vortex core radius, and eventually leads to a more com- 
plex spin structure, with the magnetization at the disk 
rim pointing down (Fig. [2]i). Further increase in the per- 
pendicular anisotropy leads to a uniform perpendicular 
magnetization, as shown in Fig. [JJd. 

The dependence of the magnetic structure of the disks 
with the value of the perpendicular anisotropy is exhib- 
ited in Fig. [3l Here we have fixed the thickness of the 
nanodot (10 nm) and obtained the magnetic phase di- 
agram in the plane of perpendicular anisotropy versus 
reduced radius of the dot derived both by micromagnetic 
simulation (red continuous line), and obtained analyti- 
cally (blue dotted line). The agreement between the two 
methods is very good. 



B. Disks: Analytical Method 



dot rim. We take into account volume (Ky) and perpen- 
dicular anisotropy {K^)^ as well as dipolar and exchange 
energy contributions. 

The energy of the magnetic states with in-plane (IP) 
and out-of-plane (OP) uniform magnetization can be 



written as 
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where is the demagnetizing factor—. 
Nz is given by^^: 
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and 



where Ji{x) is the first order Bessel function 
F2i(a,6, c, x) is the hypergeometric function. 

The vortex states can be generally described in terms 
of the magnetization Mz{p) = MsUizip)^ and it can be 
shownH that the relevant energy terms can be written as 



In order to describe the configurations of magnetic 
nanodisks we have developed a simple model for the mag- 
netic vortex state with out-of-plane magnetization at the 
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where Jq{x) are Bessel functions. For the vortex states 
we consider the following ansatz: 
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where 6 is a parameter related to the core radius^'^, c is 
related to the size of the out-of-plane magnetization at 
the rim of the dot, and ^ (0 < ^ < 1) is used to describe 
the magnetization at the rim; for an usual vortex ^ = 0. 
With the above magnetization we perform a numerical 
evaluation of the total energy with minimization of the 
adjustable parameters 6, c and g. 

This theoretical description has allowed the determina- 
tion of the phase diagram as a function of perpendicular 
anisotropy, as well as the boundaries of the region of the 
diagram where the magnetic nanodisks exhibit perpen- 
dicular magnetization at the rim (Fig. [3]). 

The vortex core radius (rc) can be defined as the value 
at which rriz = 0.5, and then Tc = (1 — 2~^/^)6, where 
b is obtained by minimization of the energy. Using the 
magnetic parameters'^ for bulk Co, (the exchange length 
is lex = 4.93 nm), we obtain the core size in the presence 
of perpendicular anisotropy, as shown in Fig. HI in good 
agreement with Fig. 1 of Garcia et al.?^. 



C. Disks: Experimental 

The samples were produced by magnetron sputter- 
ing deposition, by means of e-beam lithography on 
SiO2/Si(100) wafers. The samples presented the same 
layer structure ([Co/i/Pt2]6/Pt6) and were distinguished 
by the Co layer thickness {h = 0.6, 0.8, 1.6, 2.0 nm) in 
the stack. We have chosen these thicknesses regarding 
a perpendicular to in-plane magnetic anisotropy transi- 
tion observed when h is increased from 0.4 to 0.8 nm. 
Each sample contained arrays of 1 pm and 2 /im-diameter 
disks. For better comprehension of the results, a contin- 
uous film sample was produced along with each of the 
structured samples by placing a resist-free wafer on the 
side of the lithographed sample in the sputtering cham- 
ber. The quality of the lithography and deposition pro- 
cess has been verified by field emission gun scanning 
electron microscopy, Dektak profilometry and Rutherford 
Backscattering Spectroscopy. 
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FIG. 4: (Color online) Dependence of the magnetic vortex 
core diameter with the perpendicular anisotropy. The black 
squares are computed with micromagnetic simulation includ- 
ing Kz] the circles (red) are experimental values obtained 
by MFM (righthand scale) of Co/Pt disks and the continu- 
ous (blue) line is obtained with the analytical model. Inset: 
MFM image of 1 fim diameter Co/Pt nanodisks showing the 
vortex cores. 
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FIG. 5: (Color online) Magnetic configurations of elliptic nan- 
odots that appear in the phase diagram of Fig. (6] A) in-plane 
quasi-uniform magnetization, B) out-of-plane quasi-uniform 
magnetization, C) single vortex configuration, D) double vor- 
tex configuration. 



For the MFM measurements we used an NTEGRA 
Aura MFM scanning probe microscope (NT-MDT Co.) 
with a commercial MFM tip (NSGOl type, CoCr mag- 
netic coating, NT-MDT Co.) magnetized along the tip 
axis in the field of a permanent magnet. The MFM im- 
ages were acquired in the tapping mode at room tem- 
perature. In order to avoid instrumental artifacts in the 
determination of vortex core size from the MFM image 
we kept the lift height constant for all the measured sam- 
ples. Although MFM is not the most suitable technique 
to determine quantitatively the vortex core diameter, we 
expected to obtain the trend of the vortex diameter with 
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FIG. 6: (Color online) Phase diagram for elliptical nanodots 
without perpendicular anisotropy {Kz = 0), as a function of 
reduced height h/lex and reduced semi- axis a/lex^ drawn from 
the minimum energy computed with micromagnetic simula- 
tion. The letters correspond to regions with different spin 
configurations: A) in-plane quasi-uniform magnetization; B) 
perpendicular quasi-uniform magnetization; C) in-plane vor- 
tex; D) double in-plane vortex; E) lateral vortex. The ellipses 
in every case have semi- axes in the ratio ajh — 2. 



the thickness of the Co layer (Fig. S]). 

We have determined experimentally the trend toward 
increasing vortex core diameter in the Co/Pt multilayers, 
as the perpendicular anisotropy acting on the nanodots 
is increased. The MFM measurements made on the Co 
disks to study this effect, however, do not allow the ac- 
curate determination of the vortex core diameter. They 
only allow the observation of this increasing trend, as 
shown in Fig. HJ In the figure, the vortex core diameters 
are plotted as a function of Co layer thickness or effec- 
tive anisotropy {Keff = Ky—Kz)] in the analytical curve 
Ky = 0. Fig. HJalso shows the vortex core diameters ob- 
tained analytically and shows their agreement with the 
micromagnetic simulation results. 



III. RESULTS FOR ELLIPSES 

A. Ellipses: Micromagnetic Simulations 

Following the same methodology, we have also ob- 
tained the phase diagram of elliptical nanodots with 
Kz = 0; we have simulated ellipses where the major semi- 
axis (a) is twice the minor one (b), i.e., the ellipses in 
every case have a/b = 2. As we can see from Fig. [6l in 
this case the diagram is richer than that of the disks. The 
first observation is that the vortex state only occurs for 
dimensions that are larger than in the case of the disks, 
i.e., for approximately a > 10 lex- This is so because 
the eccentricity introduces an uniaxial shape anisotropy 



FIG. 7: Phase diagram for elliptical dots with thickness 
/i = 50 nm: diagram of perpendicular anisotropy Kz ver- 
sus reduced major semi-axis a/lex- The letters correspond to 
regions with different spin configurations: B) perpendicular 
quasi-uniform magnetization, C) in-plane vortex, E) lateral 
vortex, F) modified in-plane vortex, and G) double lateral 
vortex (see Fig. [8|). 



along the major axis, which favors a quasi-uniform in- 
plane magnetization. Another state observed is the re- 
gion in the figure which corresponds to one lateral vortex, 
that also appears in some simulations for cylinder^; to 
the best of our knowledge, this had not been observed for 
ellipses. A very interesting phase of this diagram occurs 
for a > 240 lex and h > S lex^ where two vortices appear. 
This configuration has been observed experimentally by 
several authors e.g.,^^^^. 

The phase diagram for the elliptic nanodots is also 
modified by the presence of perpendicular anisotropy. Its 
effect is illustrated in Fig. [71 which shows the phase dia- 
gram for the spin configurations obtained by micromag- 
netic simulation for ellipses with thickness of 50 nm, as a 
function of the major semi-axis a, for different values of 
the perpendicular anisotropy. This diagram is more com- 
plex than that obtained for the disks with perpendicular 
anisotropy (Fig. [3]). The letters in Figs. [6] and [71 corre- 
spond to the spin configurations: A) in-plane magneti- 
zation (Fig. [SK), B) perpendicular magnetization (Fig. 
[5p), C) in-plane vortex (Fig. [5p), D) double in-plane 
vortex (Fig. [5p), E, F, G) types of lateral vortices (Fig. 
[8]). Some lateral vortex configurations also shown in the 
elliptic nanodot cross section are illustrated in Fig. [51 
They are: E) "two-domain" out-of-plane structure with 
one lateral vortex, F) modified in-plane vortex and G) 
"three-domain" out-of-plane structure with two lateral 
vortices of opposite polarization. The lateral vortices in 
E and G occur between two domains. Note that the 
contrast color of the top view of the ellipses is in the z 
axis (perpendicular to the plane) and the cross section 
is shown with contrast in the y axis, to make the vor- 
tex structures more visible; note also that although the 
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FIG. 8: (Color online) Magnetic configurations of elliptic nan- 
odots, showing the plane of the dots and the cross sections 
with lateral vortices: E) "two-domain" configuration with sin- 
gle lateral vortex; F) modified in-plane vortex and G) "three 
domain" configuration with two lateral vortices of opposite 
polarization. Configuration E occurs in the phase diagram 
for Kz — (Fig. [6|), and configurations F and G appear 
in the phase diagram with perpendicular anisotropy, on the 
plane Kz x a. 



change their shape, for example, from circular to ellip- 
tical, or to vary their dimensions. In this work we have 
explored a different way of accomplishing this: the in- 
troduction of a perpendicular anisotropy term. We have 
observed that this leads to important modifications in 
the phase diagrams for these nanostructures, as demon- 
strated through results obtained by micromagnetic sim- 
ulation and analytical formulation. MFM measurements 
confirmed the trend of increasing vortex core diameter 
with increasing perpendicular anisotropy. 

In this work we have studied the different magnetic 
configurations of circular and elliptical nanodots, pre- 
senting them using x r phase diagrams obtained using 
micromagnetic simulation. In the case of circular nan- 
odots, a phase diagram was also obtained using an ana- 
lytical method; it agrees with the micromagnetic simula- 
tion. Measurements using the MFM technique show the 
same qualitative behavior in the dependence of the vortex 
core diameter with perpendicular magnetic anisotropy. 
The phase diagrams are also drawn for nanodots present- 
ing a perpendicular anisotropy term, and exhibit impor- 
tant differences from the = case: the region of the 
diagram corresponding to a magnetization perpendicu- 
lar to the plane increases, the region of M parallel to the 
plane is reduced, and more complex spin arrangements 
appear. 

The results presented here on the elliptical nanodots 
reveal the complexity of their magnetic behavior. In the 
range of ellipse sizes studied, several configurations ap- 
pear: in-plane and out-of-plane quasi-uniform states, one 
and two- vortex states, as well as configurations with lat- 
eral vortices. The latter structures obtained with the 
perpendicular anisotropy term are more complex and had 
not been investigated before; their detailed properties re- 
main to be studied. 

Investigations that allow the mapping of these differ- 
ent magnetic configurations are useful in designing ex- 
periments to study the basic properties of these novel 
magnetic structures, or tailoring them for technological 
applications, such as magnetic random access memories. 



ellipses in Fig. [8] are shown with the same size, they 
correspond to different semi-axes and anisotropics in the 
phase diagram of Fig. [71 

IV. CONCLUSIONS 

There are several ways of playing with the magnetic 
configurations of nanodots; the most obvious ones are to 
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